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We discuss the cosmological reconstruction in modified Gauss-Bonnet (GB) gravity. It 
is demonstrated that the modified GB gravity may describe the most interesting features 
of late-time cosmology. We derive explicit form of effective phantom cosmological models 
ending by the finite-time future singularity (Big Rip) and without singularities in the future 
(Little Rip). 



Recent observational data indicate that our Universe is currently in accelerated phase pQ. Since 
the discovery of the acceleration, a large number of possible mechanisms have been proposed to 
explain the origin of the dark energy. General Relativity in its standard form can not explain the 
accelerated expansion without extra terms or components, which have been gathered under the name 
of dark energy. There are various proposals to construct an acceptable dark energy model, such 
as: a scalar field with quintessence-like or phantom-like behavior [2], spinor, (non-)abelian vector 
theory, cosmological constant, fluid with a complicated equation of state, and higher dimensions 
[3]. As a result, we have a number of competing scenarios which to describe late-time acceleration. 
Working in framework of one of the above proposals, we are forced to introduce the following extra 
cosmological components: inflaton, a dark component and dark matter. 

We concentrate on the modified gravity [U |5], which represents a classical generalization of 
general relativity (modifications of the Hilbert-Einstein action by introducing different functions of 
the Ricci scalar [5j |7] or Gauss-Bonnet invariant [8, 9j), should consistently describe the early-time 
inflation and late-time acceleration, without the introduction of any other dark component. 

In the present work we will focus on the Gauss-Bonnet gravity, where the gravitational action 
includes functions of the Gauss-Bonnet invariant 



I. 



INTRODUCTION 



G = R 2 - iR^R^ + R^xaR^ ■ 
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It is interesting, that the theory of this type is closely related to higher order string curvature 
corrections |10j . The starting action is 

R 1 



d XyJ—g 



2k 2 2 



d^t-VW-tiWG 



Here G is the Gauss-Bonnet invariant ([I]) and the scalar field <j> is canonical in 
One may redefine the scalar field (j) by <fi = e<p. The action takes the following form 



S 



d x^J—g 



Here 



If a proper limit e — > exists, the action ^ reduces to 

R 



S = Id xy— ~g 



2k 



Then <p is an auxiliary field. By variation of <p, we find 



, V(<p) - £i((p)G 



o = v'(v)-£MG , 



which may be solved with respect to ip as 



p = $(G) 



Substituting Q into the action ([5|, we get for the i ? (G ! )-gravity 

R 



d X\f— g 



2k 



, F(G) 



, F(G)^vmG))-^(<i>(G))G 



(2) 



(3) 



(4) 



(5) 



(6) 



(7) 



(8) 



On the other hand, the action ?? is associated with non-local models Let's start with the 
next action 

R K 2 



S = Id X\f^g 



GU~ X G + L r 



2k 2 2a 



(9) 



Here C m is the matter Lagrangian. By introducing the scalar field <j), one may rewrite the action 
([9]) in a local form: 

R a 



S 



2k 2 2K 2&iM 



+ (t>G + £, 



(10) 
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In fact, from the ^-equation we find = — ^□ _1 G. By substituting this expression into (10), we 
reobtain (|9]). 

Thus, by studying F(G) - gravity, we are also exploring several other models. Note that in some 
cases, F(G) or even F(R,G) gravity maybe transformed to F(R) gravity |12] , 

The five-year Wilkinson Microwave Anisotropy Probe (WMAP) data [T] give the bounds to the 
value of the equation of state (EoS) parameter wde, which is the ratio of the pressure of the dark 
energy to the energy density of it, in the range of -1.11 < wde < -0.86. This could be consistent 
if the dark energy is a cosmological constant with wde = -1 and therefore our universe seems to 
approach asymptotically de Sitter universe. 

Although the accelerating expansions seem to be de Sitter type, the possibility that the current 
acceleration could be quintessence type, in which wde > -1; or phantom type, in which wde < -1) 
is not completely excluded. 

It is often assumed that the early universe started from the singular point often called Big Bang. 
However, if the current (or future) universe enters the quintessence/phantom stage, it may evolve 
to the finite-time future singularity depending on the specific model under consideration and the 
value of the effective EoS parameter (Big Rip) [6j. An elegant solution to this problem has been 
recently proposed |14j under the form of the so-called Little Rip cosmology which appears to be a 
realistic alternative to the ACDM model. 



II. [R + f(G)} GRAVITY 
Consider the following action (see Ref. [9]): 



S = I (fx^g ( R + f{G) + L m ) (.1.1) 



where k 2 = 8itGn, Gn being the Newton constant, and the Gauss-Bonnet invariant is defined 
by([T]). By varying the action (11) over g pv , are obtained: 



= (-R^ + \g^R\ + T>» + \ffiG) - 2f G KRT + 4f G R»R va - 2f G R^R^ T 

- Af G R^R aP + 2 (V»V»f G ) R - 2g^(V 2 f G )R - A(V p V^f G )R^ 

- 4(V P V7 G )^ + 4(V 2 f G )R^ + Ag^(V p V a f G )R^ - 4(V a V pf G )R^ , 



where f G = f{G) and f GG = f"{G) (' = d/dG). 



4 



We consider a spatially-flat universe 

3 

ds 2 = -dt 2 + a(t) 2 J2(.dx 1 ) 2 ,. (12) 



Here a(t) is the scale factor at cosmological time t. It is easy to obtain the FRW equations 
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= -^H z + Gf G - /(G) - 24GH foG + Pn 



= 8H 2 f G + 16H(H + H 2 )f G + ^(2H + 3H 2 ) + / - Gf G +p m . (13) 



1 

K 

Here H is the Hubble rate {H = a/a,'= d/dt), p m is matter energy density: 



p m + 3H(1 + w)p m = , (14) 
while G, its derivative and R can be defined as functions of the Hubble parameter as 

G = 2A{HH 2 + H i ), R = 6(H + 2H 2 ), 
G = 24H(HH + 2H 2 + AH 2 ). 



It is easy to see that the second equation in ( 13 ) is a differential consequence on the solution 



The first equation (13) can thus be expressed as follows 



= ^H 2 + 2A{HH 2 + H 4 )f G - /(G) - 2A 2 H A {HH + 2H 2 + 4H 2 )f GG + p m , (15) 

By selecting certain of the Hubble parameter H we obtain a differential equation for the function 
/. This is the method of reconstruction (this method has been implemented in Ref. [13J for f(R) 
gravity and Ref. [9j for f(G) gravity). 

III. f(R, G) GRAVITY 

Let us now consider a more general model for a kind of modified Gauss-Bonnet gravity. This 
can be described by the following action 

S = d A Xyf^g 



(16) 



^F(R,G) + L % 

Varying over g pv the gravity field equations are obtained |9], 

= T» v + ^g^F(G) - 2F RtG RR" u + 4F G R p l R up 
-2F G R?P° T R v p(JT - AF G R ppau R pa + 2{V P V U F G )R - 2g pv {V 2 F G )R 
-4(V p V p F G )R up - 4(V p VF G )R pp + A(V 2 F G )R pu + 4g pu (V p V a F G )R pa 

-A(V p V a F G )RV pua - F R R pu + V"V U F R - g pu V 2 F R . (17) 
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In the case of a flat FRW Universe, described by the metric (??), the first FRW equation yields 

= ^{GFg -F- 2AH 3 F Gt ) + 3(H + H 2 )F R - 3HF Rt + k 2 p m . (18) 
For simplicity, we consider the following subfamily of functions 

F(R,G) = f 1 (G) + f 2 (R) . (19) 



Correspondingly, the Friedmann equation (18) can be split into two equations, as 



= -2AH 3 Gf 1GG + G/ig - /i , 
= -3HRf 2RR + 3{H + H 2 )f 2R - l -h + K 2 p m (20) 



IV. LITTLE RIP MODEL 

In Ref. |15| is proposed Little Rip model on the basis of a viscous fluid. Using the obtained in 
this study form the scale factor, we construct a model of the GB gravity. 
Consider the scale factor as |15 



a{t) = expa [e 131 - I J (21) 
In this case the Hubble parameter has the form 

H{t) = a/3e pt (22) 
Then the derivative of the Hubble parameter is proportional to itself, and 

G = 2AH 3 (j-o^j. (23) 
FRW equation for the R + f(G) theory gravity takes the form 



n _ 3 „ 2 , , 9(3 2 H + 19(3H 2 + 4H 3 df (3H 2 d 2 f 

U " k* H /+ (3/3 + 4F)2 dH 3(3 + 4H dH* + Pm - [M) 
Choose the energy density in the form p m = const we obtain the following solution 
18f3 3 H + 3(3H 3 + 3H 4 + 2(3 2 (6H 2 - K 2 Pm ) H 3 (P + H)d J 

HC 2 (-PeT (3/3 2 + 2(3H + H 2 ) + H 2 (f3 + F)ExpIntegralEi 



+ 



2/33(1 + p) 



I. 




3xl0 18 - 



FIG. 1: Plot of f(t), a) (t,0,l), b) (t,0,10) (Green line - C 2 = 0, C x = 10, red line - C 2 = 0, d = 0, blue 
line - C 2 = 10, C\ = 10, black line - C 2 = 10, C x = 0). 



Here ExpIntegralEi[z] gives the exponential integral function where Ei(z) = 

This solution can be illustrated by the graph (we assume that a=l, (3 = 1k = 1, p m = 0): 
Figures show that at the initial time all the solutions have the same behavior (dominated by the 



first term in (25)), then begins to play the role of the second term (C\ ^ 0), and finally turned the 
third term [C 2 7^ 0), this is clearly seen in the figure 1. 





FIG. 2: Plot of f(t), (t,0,2.2) (Red line 
C 2 = 1, Cx = 1, blue line - f(t) = 4e e *). 



FIG. 3: Plot of /(G), (G,0,200) (Red line 
- C 2 = -10, C x = 100, blue line - f{G) = 
8G^e Gl/i ). 



We can consider the dependence of the F(G). Obviously not possible to do this, but for large H, 



we can assume that G = 2AH 4 . From (25) one can see that there are 3 terms: the first dominates 
at small G, then by 30 < G < 100 included the last term (C 2 + and f(G) ~ G 1 / 8 e Gl/ \ figure 3), 
then included the second term (C\ 7^ and f(G) ~ G, figure 4a) and then G— > 00 (figure 4b) 
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again dominated by the last term. 




2000 4000 6000 8000 10000 Z0OO0 40000 60000 80000 100000 



FIG. 4: Plot of /(G) a) (G,0, 10000), b) (G,0,100000) (Red line - G 2 = -10, C x = 100, blue line - /(G) = 
1.8G). 

This behavior is due to the fact that the last term consists of two parts with different signs. For 
small values of G terms rather different, and with increasing G, they become the same order and 
the difference becomes smaller than other terms. When G tends to infinity, the difference between 
them begins to grow faster than the other terms. 

By selecting different values of the constants C\ and C2 can be different depending / from G. 
For example, we construct the graph for the following f(G) values of the constants of integration: 
C 2 = 1, -20 < Ci < 20 (figure 5). 




FIG. 5: Plot of /(G), (G,0,70) (G 2 = 1, -20 < d < 20). 
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Consider now the F(G, R) - gravity. Let us solve the equations (20) 
h(H) ' 



(4a(2 + a(2 + (-5 + a)a))/3 4 ) 
(^e^H (/3 2 + 4/3# - H 2 ) C 2 + H (2/3 3 + 2fi 2 H - 5(3H 2 + H 3 ) (^d + C 2 ExpIntegralEi j J 



(26) 



h(H) 



2a 3 (l + a)/3 4 



(27) 



-a 2 f3e^H (3/3 2 + 2/3H + H 2 ) C[2] + H 3 ((3 + H) I 2C[1] + a 2 C[2]ExpIntegralEi 
Using the 

# = ~ (-3/3 + V9/3 2 + 12fl) • 
We obtain the dependence of the function f 2 from i? 



if 

1 



(28) 



/ 2 (i?) = (-3/3+^9/3 2 + 12i?) ( -6(3e~ 30W ^^^ (2R + 3/3 f/3 - 3x/9/3 2 + 12i?)) C 2 + 



117 P 3 - 69/3i? + 57/3 2 v / 9 ( S 2 + 12i? + iV9/3 2 + 12i2) 

3/3 + ^9/3 2 + 12i2" 



(29) 



4Ci + C 2 ExpIntegralEi 



12/3 



We can consider the dependence of the fi(G). Obviously not possible to do this, but for large 
H, we can assume that G = 2AH 4 . Then 



G V4 



fi(G) = -6 a \ 2 p e ^m ( 12V3/3 2 + 4 e 1 / 4 ^ 1 / 4 + V2VGJ G x ^C 2 + 

G l/4 



+ ( 6\/2/3G 3/4 + 6 3/4 G) 2Ci + al 2 C 2 ExpIntegralEi 



2 3 / 4 3V4/3 



(30) 



V. POWER LAW SOLUTIONS 



Consider a power-law solutions. Choose the Hubble parameter as 



Hit) 



a 

t„-t' 



(31) 
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Then the scale factor is given by 

ait) = jj-ljp (32) 
This scale factor can meet the various phases of expansion, which are characterized by different 



values of the equation of state (EoS) parameter wde- To easily find the value wde f° r (32). The 
equations FRW have the form: 

P=^2 H ^ P = -^{ 2 H + 3H"). (33) 



Comparing with the (13), we can write formally the same equations by removing all unnecessary to 



the effective pressure and energy density;: 



So 



Peff = Gf G - f(G) - 2AGH 3 f GG + Pm 

p ef f = 8H 2 f G + 16H(H + H 2 )f G + f-Gf G + Pm . (34) 



p eff 1 2H + 3H 2 

W DE = -^ = -- 772 • ( 35 ) 



For (|31J have w DE = -1/3(2/q + 3). 

Consider the different values of wde- If Wde = -1 this could be consistent if the dark energy is 
a cosmological constant and therefore our universe seems to approach to asymptotically de Sitter 
universe. This value is not possible for the chosen scale factor. This value is the min in the 
permissible region, but the possibility that the current acceleration could be quintessence type, in 
which wde > -1> or phantom type, in which wde < -1> is not completely excluded. For quintessence 
type, the value of a is negative. For phantom type, the value of a is positive. 



It is easy to write of the equation (15) for this cases as an equation for f(t) 



= -4(1 + a) (3a 2 - K 2 p (t s - t) 2 ) + 

+ K 2 {t s - t) 2 (-4(1 + a)f(t) - (t s - t) (-(6 + a)f(t) + (t, - t)f"{t))) (36) 



where pq = const. 



(37) 
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We now express t through G 



23/4 3 l/4( a 3 +a 4)l/4 

ts-t- ^yl • ^») 



And we obtain the following form of the function f(G) 



= _ V iV" 3 (l + Q0 g _ ._ 2a2 _a (1+a)3 i ( _ 1 _ a)G i± ?s ,3 + , |(-l-a) ^ (39) 

(—1 + aja^r 

It is seen that if a = 3 the equation is considerably simplified 

3^ 



y/2. 



K 

Now consider the space filled with, dust p m — ^-jj . 



po. (40) 



_ v fV« 3 (l + a)G 2 2+g f 3 3 °/4 (1 + a)G -3a/2 ( q3(1 + a )G) 3Q / 4 po 
/( )_ (-l + a)oK 2 + (4 + 19a + 12a 2 )C 3 

_ 2 -i( 1 +«) 3 |(-i-«) e -«- G ^ (a 3 (l + a)G) * (_1 ~ a) C x (41) 

It is seen that if a = 3 the equation is considerably simplified 

f(r\ 644972544 2*/% 

/l J " y^ 2 i69C 3 G9/4 ' 



Consider now the F(R, G) - gravity (20). Equations take the form 

= 4(1 + a)h - (t. - t) ((6 + a)fi - (t s - t)ft), (43) 



= (2 + 4a) f 2 - (t„ - t) ((4 + a)f 2 - (t s - t)%). (44) 

The solution of these equations have the form 

{t _ ts) s- aCl + C 2 2 3 / 4 3 V 4 (a 3 G + a 4 G) 1/4 
/l = ' ts " t = ^ (45) 



/ 2 = (t-t 5 )K- 3 — Vi + (-io +Q)a ) f Cl + i t-t s )V^(-^c 2 ),t s -t= " MVa ^ +2a) . (46) 

^ ' v R 

Or finally get 

/i(G) = CiG (1+Q)/4 , (47) 
/ 2 (fl) = J2 K 3+ * f Vi+(-i°+«)«) + jj-IVi+C-io+^-Ca) . (48) 
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VI. DE SITTER SOLUTIONS 

De Sitter solutions are described by an exponential expansion of the Universe, where the Hubble 
parameter and the scale factor are given by 

H(t) = H ^ a(t) = e Hot , (49) 

where Hq is a constant. This kind of solutions are very important, as the observations suggest that 
the expansion of our Universe behaves approximately as de Sitter. It has been shown in Ref. |16| 
that de Sitter points are critical points in f(R) gravity. It is straightforward to see that this is also 
the case in R + /(G) gravity. 

= -^Hl + Gof G (G ) - /(Go) • (50) 

Here, Go = 2AHq and we have ignored the contribution of matter. Then, we have reduced the 
differential equation to an algebraic equation that can be resolved by specifying a function /(G). 



Consider the equation (50) as a differential equation for G: 



3\/G 
2 2^ 



° = -\/o^2 +G/ G (G)-/(G) . (51) 



In this case, the function F{G) takes the form 



/(G) = (52) 



We see that the function of the form /(G) = cyG describes the power solutions and de Sitter 
solution (c = The general solution has the form 

a(t) = (t-at-Ci)^C 2 , (53) 

here C\ and G2 -constants of integration. The case when a = 1 corresponds to the de Sitter solution. 
For F(G, R) gravity similarly, we have the following equations 

= ~\ (G F G (G ) - F(G , Ro)) + 3H$F R (R ) . (54) 

If we consider the case F(G, R) = /i(G) + /2(-R) is not difficult to find a possible solution F(R, G) = 
R 2 . 
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VII. LITTLE RIP IN GAUSS-BONNET GRAVITY WITH DILATON 

We can find a Little Rip solution for Gauss-Bonnet gravity with dilaton for spaces of dimension 
D + 4. Consider space oiD=p + q + l dimensions with two maximally symmetric subspaces: 
p-dimensional and g-dimensional. We write the Lagrangian as 

L = ^2~ \%V9 V V - V{<p) + e(<p)G. (55) 

Here £(<f), V((p) are functions of dilaton ip. We choose them as 

e(<p) = /3e" w , V(tp) = ae™. (56) 



Variating the action with Lagrangian ( 55 ) we get field equations [3] . One of the possible solutions 
of these equations will be 

a(t) = a exp|^ e ^/ 2 |, (57) 

6(t) = 6 exp| 2 ^°e^/ 2 | (58) 

where c*o, &o and cq are arbitrary positive constants, 1*1,1*2, 1*3 are constants satisfying the following 
equations 

3^i 2 + ^qiU2 2 + 3gi*iit 2 + ^ii3 2 ~\ a + ^{36gini 2 u 2 2 + q 3 u 2 4 + 

+ \2u\u 2 (2qu 2 + q 2 u 2 2 )} - 2(3ju 3 {18qu 1 2 u 2 + 9qiUiu 2 2 + 6gi*i 3 + q 2 u 2 3 } = 0; 
2 l t„ , n„„. 2 / , \ „„, ...la," 1 , , . 2 



-3ui 2 - -(g + l)gi* 2 - 71*3 ^1 + -i* 2 J - 2qu x u 2 + -« 3 + - - -^{-4c7(3*7 + l)7i* 3 i*iu 2 
- 28<nii 3 ui 2 -U2 - 2(q + 2)<?i7u 3 u 2 3 - ^ 2 u 2 u 2 - 2q 1 j 2 u 3 2 u 2 2 + Aq(5q - 3)ni 2 u 2 2 + 



+ 8*7*71 1*11*2 + 16*7i*i 3 i*2 + {g + 1) 1721*2 — 167-t*3iii 3 — 8*77 113 1*11*2} = 0; 

In 2 1 ~ 1 ^ lo 

- 2 91^2 -61*1 — 7113112 - -71*31*1 - 6{q-l)u\u 2 + -u 3 

- ^(3{-2(q + l)(q- l) 2 ^u 3 u 2 3 - 6(q - l)(3q - 2)^u 3 uiu 2 2 - 60(q - l)ju 3 ui 2 u 2 - 

- 60 7 i*3i*i 3 - 2(q - 1) 2 7 W - 12 7 2 u 3 2 i*i 2 + q 3 u 2 A + 2A(q - l)(2q - 3) Ul 2 u 2 2 - 



+ I2(q - l){q - l) 2 niii 2 3 + 72{q - l)i*i 3 u 2 + 24ni 4 - \2{q - l)7 2 u 3 2 ixi« 2 } = 0; 
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-7U 3 2 - 2(3m + qu 2 )u 3 + /3<y{12q(4q - 2)u 2 u 2 2 + {q + l)q 2 u 2 4 + 72qui 3 u 2 + 12qqiuiu 2 3 + 

+ 18g i i7ii3iiiM2 2 + 3Gq~/U3Ui 2 u 2 + I27U3W1 3 + 24ui 4 + 2(727ii 3 «2 3 } + 07 = 0. 
Therefore it is necessary to find solutions of this equations satisfied conditions 

Ui > 0, u 2 < 0, U3 > 0. 

This is easily done, for example 

q = l : a = 1, /3 = 1, 7 = 1, u 3 = 0.383, m = 0.378, u 2 = -1.32. 

g, = 2, a = 0.01, = 100, 7 = 0.01, u 3 = 5.09 • 10~ 5 , u x = 0.0175, u 2 = -0.0803. 
q = 3 ] a = 0.001, (3 = 100, 7 = 0.001, n 3 = 2.33 • 10~ 4 , m = 0.0306, u 2 = -0.0788. 

VIII. CONCLUSIONS 

In this paper we construct a series of cosmological models which describe the accelerated expan- 
sion of the universe. The models in framework of the modified Gauss-Bonnet gravity are constructed 
by the reconstruction method. Constructed models are of the effective phantom type and do not 
lead to a singularity in the future (Little rip models). 

An explicit form of the action leading to such a model within the framework of gravity R + f(G) 
and f 2 {R) + fi(G) is obtained. The solution contains special function, which can be removed by 
choosing the constant of integration. The resulting action consists of three terms, each of which 
plays a role at the definite stage of the evolution of the Universe. Since the last term represents the 
difference between the functions of the same order, and this order of these functions is maximal, 
it has double the lead - with the average values and tends to infinity. In addition we construct a 
model describing the Little Rip in the multivariate Gauss-Bonnet theory with dilaton, when visible, 
and additional spaces are maximally symmetric. 

Little Rip models provide an evolution for the universe intermediate between asymptotic de 
Sitter expansion and models with a big rip singularity. Therefore, the method of reconstruction 



models were built describing the de Sitter space and space with the metric 32 which refers to the 
phantom type and has a singularity of the future. 

For all constructed models wde < — 1 (except deStter model) and therefore they belong to the 
phantom type. For models of the effective phantom energy conditions are always broken. Such 
behavior is typical for the Little Rip models |17] 
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For Little Rip models with Hubble parameter (22), we obtain 

2af3 2 e^ 



P + P 



K 2 



We can continue to explore these theories and consider stability of solutions. However, as shown 
by work [18J Little Rip solutions are stable and can be expected and received from our paper 
solutions will have the same property. 
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